Spectrum of extended systems from Reduced Density Matrix Functional Theory 
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We present a method for calculating the spectrum of extended solids within reduced density matrix 
functional theory. An application of this method to the strongly correlated transition metal oxide 
series demonstrates that (i) an insulating state is found in the absence of magnetic order and, in 
addition, (ii) the interplay between the change transfer and Mott-Hubbard correlation is correctly 
described. In this respect we find that while NiO has a strong charge transfer character to the 
electronic gap, with substantial hybridization between ti g and oxygen-p states in the lower Hubbard 
band, for MnO this is almost entirely absent. As a validation of our method we also calculate the 
spectra for a variety of weakly correlated materials, finding good agreement with experiment and 
other techniques. 
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A derivate of the ground-state density functional the- 
ory (DFT) calculations are the Kohn-Sham (KS) eigen- 
values, which lead to a non-interacting spectrum. Even 
though the KS equations represent an auxiliary non- 
interacting system whose states and eigenvalues may be 
quite different from the true quasi-particle system, em- 
pirical evidence shows that in many cases this single 
particle KS spectrum is in agreement with the x-ray 
photo-emission Spectroscopy (XPS) and Bremsstrahlung 
isochromat spectroscopy (BIS) experiments. However, 
for strongly correlated materials, this KS spectrum is in 
fundamental disagreement with experimental reality. In 
the absence of spin-ordering all modern exchange cor- 
relation (xc) functionals within DFT fail to predict an 
insulating ground-state for transition metal mono-oxides 
(TMOs), the prototypical Mott insulators. On the other 
hand, it is well known experimentally that these materi- 
als are insulating in nature even at elevated temperatures 
(much above the Neel temperature) [l|, |2| . This indicates 
that magnetic order is not the driving mechanism for the 
gap, but is instead a co-occurring phenomenon. In fact 
not only DFT, but most modern many-body techniques 
also fail to capture the insulating behavior in TMOs with- 
out explicit long range spin ordering. 

In this regard reduced density matrix functional theory 
(RDMFT) has proved to be valuable in that it not only 
improves upon the KS band gaps for insulators in general, 
but also predicts TMOs as insulators, even in the absence 
of long range spin-order [3j. This clearly points towards 
its ability to capture physics well beyond the reach of 
most modern day ground-state methods. Despite this 
success the effectiveness of RDMFT as a ground state 
theory is seriously hampered by the absence of a tech- 
nique for the determination of spectral information. 

In this work, we present a method for calculating the 
spectrum of extended systems within the framework of 
RDMFT. Using this approach the spectra of a set of 
TMOs is determined. In all cases we find good agree- 



ment with experimental data as well as with the results 
of GW and dynamical mean field theory (DMFT) calcu- 
lations. As an additional test we determine the spectra 
for a set of systems without strong correlations; the met- 
als Cu and V and the band insulator BN. The spectra 
thus found show good agreement with experimental data 
as well as with data from other techniques, indicating a 
general applicability of the method proposed here. 

Within RDMFT the one-body reduced density matrix 
(1-RDM) is the basic variable 0, i] 

7 (r,r') = N J dr 2 . . . dr N $*(r', r 2 . . . r N )$(r, r 2 ...r N ), 

(1) 

where $ denotes the many-body wave function. Di- 
agonalization of this matrix produces a set of natural 
orbitals[J], (f'j'k-i an d occupation numbers, rijk, leading to 
the spectral representation 
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where the necessary and sufficient conditions for ensem- 
ble 7V-representability of 7 '(| require < njk < 1 for 
all j, k, and X)j k n i k = Here j represents the band 
index and k the crystal momentum. 

In terms of 7, the total ground state energy [5| of the 
interacting system is (atomic units are used throughout) 



£M=-~ [ lmxVh(r,r')d 3 r' + J p(r)V c ^(r)d 3 



e^d^r' + EM 
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where p(r) = 7(1-, r), V CKt is a given external potential, 
and E xc we call the xc energy functional. In principle, 
Gilbert's [B| generalization of the Hohenberg-Kohn the- 
orem to the 1-RDM guarantees the existence of a func- 
tional E[y] whose minimum yields the exact 7 and the ex- 
act ground-state energy of systems characterized by the 
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external potential V BX t(r). In practice, however, the cor- 
relation energy is an unknown functional of the 1-RDM 
and needs to be approximated. Although there are sev- 
eral known approximations for the xc energy functional, 
the most promising for extended systems is the power 
functional [3j where the xc energy reads 



E xc [l] = #xc[Wik},{roik}] 
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In this functional we set a = 0.565, so that the corre- 
lation energy of the homogeneous electron gas (HEG) 
is well reproduced in the range of r s appropriate for 
solidsQ- The method then becomes truly ab-intio in 
the sense that the theory does not depend on parameters 
that are adjusted to the system to be calculated. 

The main aim of the present work is to determine the 
ionization potentials and electron affinities of extended 
systems within the framework of RDMFT. To achieve 
this the total energy functional is first minimized with 
respect to {</>i k } and {n^} for the -/V-particle system to 
obtain the ground state. The electron removal energies 
and electron addition energies arc then calculated as 



dE[{</>i k }, {n ik }} 



«ik=l/2 



(5) 



Numerical tests show that the dependence of the to- 
tal energy upon varying a single occupation number rij k 
is essentially linear, see Fig. [TJ and thus Eq. where 
the derivative is calculated at half filling, describes the 
change in energy when a particle with momentum k is 
either added or removed from the system with all other 
degrees of freedom held fixed. For Q k > n (/j, being the 
chemical potential) this derivative then corresponds to 
the electron addition energy (En+i—En) and for e^ < fi 
to the electron removal energy (En — Ejsr-i). 

We now deploy this scheme for the determination of 
the spectra of several extended solids; the strongly cor- 
related Mott insulators NiO, CoO, FeO and MnO; the 
band insulator BN; the transition metal V, and the Nobel 
metal Cu. All calculations are performed using the full- 
potential linearized augmented plane wave code Elk[8|], 
with practical details of the calculations following the 
scheme described in Ref. 0. 

Presented in Fig. [2] are the spectra generated via Eq.[5] 
for the Mott insulators under consideration. Also shown 
are GW data taken from Refs.l9l.Hol and dynamical mean 
field theory (DMFT) results that for NiO taken from 
Refs. [H Gj and for MnO from Ref. 13. For details of 
these calculations we refer the reader to the afore men- 
tioned works, however we note that in all cases the cal- 
culations are spin polarized. The experimental data are 




FIG. 1: Change in total energy upon changing a single oc- 
cupation number n ik , with all others held fixed. Inset shows 
the zoom of the same for small values of n; k - 
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FIG. 2: (Color online) Density of states (DOS) for the tran- 
sition metal oxides NiO, FeO, CoO, and MnO. Shown are 
XPS and BIS spectra, in addition to calculations using the 
GW, DMFT, and RDMFT methods. The GW and DMFT 
results are from spin-polarized calculations, and are vertically 
shifted for clarity, while the RDMFT calculations are spin- 
unpolarized. 



from Refs. Q [B, E3 for MnO, CoO, NiO and FeO 
respectively. 

It is immediately apparent from Fig. [2] that RDMFT 
captures the essence of Mott-Hubbard physics: all the 
TMOs considered are insulating in the absence of spin 
order. Furthermore, examination of the spectra for NiO 
and CoO reveals an excellent agreement between the 
RDMFT peak positions and the corresponding XPS and 
BIS data. In fact, not only are the peak positions well 
reproduced, but also their relative weights. Turning to 
MnO, one notes a similarly good agreement between 
experiment and the RDMFT results for peak positions 
with, however, somewhat worse agreement regarding the 
relative weights. 

In the case of FeO it must be recalled that Fe segrega- 
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FIG. 3: (Color online) Partial density of states for TMOs NiO, 
FeO, CoO, and MnO. In all cases the gap is between states of 
e g and t2 g character. Oxygen-p states are also presented to 
highlight the hybridization effects. 



tion, unavoidable in this compound, precludes the exper- 
imental realization of pure FeO samples. For this reason 
the only existing experimental data is rather old, and the 
presumably substantially broadened data presents no dis- 
tinct features that may be used for comparison. Turning 
to a comparison of the RDMFT spectra with the corre- 
sponding GW and dynamical mean field theory (DMFT) 
results, one notes that while for NiO all three methods 
are in close agreement, for MnO the peak positions are 
more closely matched between the GW and RDMFT cal- 
culations. 

As is well known, while the insulating state of TMOs 
is driven by a charge localization due to strong Coulomb 
repulsion (Mott-Hubbard correlation^ an important aux- 
iliary mechanism is charge transfer [18j due to hybridiza- 
tion between ligand and transition metal (TM) states. 
Amongst the TMO series this latter mechanism is gener- 
ally believed to play an important role in the case of NiO, 
but to be of decreasing importance as atomic number is 
lowered, with the insulating state of MnO thought to be 
driven entirely by Mott-Hubbard correlation. Clearly, an 
outstanding challenge for any ab-initio theory is to cap- 
ture both these aspects of TMO physics. 

In Fig. [3] we present the site and angular momentum 
projected DOS for the TMOs considered in this work. 
The electronic gap, as expected, always occurs between 
lower and upper Hubbard bands dominated by states of 
t2 g and e g character respectively. However, while for NiO 
one finds a significant component of oxygen-p states in 
the lower Hubbard band, for the other TMOs this hy- 
bridization between oxygen-p and TM-d states is much 
reduced, and is almost absent in the case of MnO, indi- 
cating that for this material the insulating state is driven 
mostly by Mott-Hubbard correlations. Thus we find that 
not only does RDMFT capture the Mott-Hubbard charge 
localization in the TMO series, but in addition also the 
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FIG. 4: (Color online) Difference between LSDA charge 
density and charge density calculated using LSDA+t/ and 
RDMFT, (p(r) — Plsda{*)) for CoO. Positive values indicate 
localization of charge as compared to LSDA. 



subtle interplay with charge transfer effects. 

A confirmation of this charge localization may be seen 
in the charge density difference p(r) — Plsda(t), shown in 
Fig.Hfor RDMFT and LSDA+J7 calculations of CoO. A 
comparison with LSDA+f/ is instructive as this method 
(with an appropriate choice of U) is able to accurately re- 
produce the insulating gaps of the TMO series. As CoO 
has an odd number of electrons per unit cell, spin polar- 
ization is required for the LSDA+t/ calculation to find an 
insulating ground state. However, despite the fact that 
both the RDMFT and LSDA+J7 calculations yield simi- 
lar gaps, a significant charge localization is seen only for 
the RDMFT calculation. Interestingly, one observes an 
almost spherical charge accumulation at the oxygen site, 
a result in agreement with experiment [19(, but different 
from that found in the corresponding LSDA+C/ result. 

The actual values of the insulating gaps that may be 
extracted from Fig.[2]are 4.0 eV (4.3 eV), 2.5 eV (3.4 eV), 
2.4 eV (3.4 eV), and 2.0 eV (3.9 eV), for NiO, CoO, FeO, 
and MnO respectively. The corresponding experimental 
gap is given in parenthesis. One notes that the agreement 
is closest for the NiO which has the lowest magnetic mo- 
ment of 1.7 hb amongst the TMOs considered here, with 
the worst agreement for MnO which has the largest mo- 
ment of 4.7 [Mb- As the RDMFT calculations presented 
here are non-magnetic (i.e., spin degenerate) the trend is 
natural, and indicates that for the large moment TMOs 
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the co-occurring magnetic order does contribute signifi- 
cantly to the magnitude of the gap. 



Vanadium 



counted in this approach. It should be stressed that all 
materials are calculated using the same RDMFT func- 
tional that is designed to to reproduce the correlation 
energy of the electron gas at typical metallic densities, 
i.e. the approach is universally applicable to solids. 
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FIG. 5: (Color online) DOS for (a) Vanadium, (b) Copper 
and (c) cubic-BN in comparison with XPS and BIS spectra. 
Note that the LDA data is scissors corrected for BN. The 
experimental data for BN is from Ref. [2fJ and for Cu and V 
from Refs. H HI The GW data for Cu and V is from Ref. 

m 

An essential validation of the RDMFT approach 
demonstrated here consists of showing that the physics 
of strongly correlated TMOs is not captured at the ex- 
pense of poorly treating weakly correlated systems. In 
Fig. [5] are shown the DOS for the band insulator BN, as 
well as the metals V and Cu. All these are systems with 
weak correlation and, reassuringly, we see that V and Cu 
are indeed found to be metallic in RDMFT, while the 
spectrum of BN is found to be in good agreement with a 
scissors corrected LDA calculation and experiment. The 
gap for BN in RDMFT is found to be 6.9 eV, close to 
the experimental value of 6.39 eV. The metallic DOS in 
RDMFT also shows broad agreement with experiment 
and other calculations. 

To conclude we have presented a method for the calcu- 
lation of spectra within the framework of RDMFT based 
on the derivative of the total energy with respect to oc- 
cupation number at half filling. Owing to the linear be- 
haviour of the total energy upon varying a single oc- 
cupation number, this derivative represents the change 
in energy upon addition or removal of a particle. We 
have shown that the spectral information obtained in 
this way gives a detailed account of the strongly corre- 
lated nature of the TMOs, including the balance between 
Mott-Hubbard correlation and charge-transfer character 
in these materials. Furthermore, both metals as well as 
weakly correlated band insulators are found to be well ac- 
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